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Abstract 



Exact numerical diagonalization of the Bohr Hamiltonian by SU(1, 1) x SO(5) methods is used to obtain detailed quantitative 

ON predictions for single-phonon and multi-phonon excitations in well-deformed rotor nuclei. Dynamical y deformation is found to 

significantly influence the predictions through its coupling to the rotational motion. Basic signatures for the onset of rigid triaxial 

~\ deformation are obtained. 
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The Bohr collective Hamiltonian has served as a conceptual 
benchmark for the interpretation of quadrupole collective dy- 
namics in nuclei for several decades [1,2]. A tractable scheme 
for numerical diagonalization of the Bohr Hamiltonian, the al- 
gebraic collective model (ACM) [3-7], has recently been pro- 
posed, based on SU(1, 1) x SO(5) algebraic methods. The need 
for such an approach arises since the conventional approach to 
numerical diagonalization of the Bohr Hamiltonian, in a five- 
dimensional oscillator basis [8-10], is slowly convergent and 
requires a large number of basis states to describe a general de- 
formed rotor- vibrator nucleus. Consequently, it has been nec- 
essary to apply varying degrees of approximation in addressing 
the dynamics of transitional and deformed nuclei, as in the clas- 
sic rotation-vibration model [11] and rigid triaxial rotor [12] 
treatments of the Bohr Hamiltonian, or in more recent studies 
of critical phenomena [13-16]. 

The ACM scheme, in conjunction with recent progress in 
construction of the relevant SO(5) D SO(3) Clebsch-Gordan 
coefficients [7], now permits the diagonalization of the Bohr 
Hamiltonian for potentials of essentially arbitrary stiffness, as 
considered in this letter. The Bohr Hamiltonian can conse- 
quently be applied, without approximation, to the full range of 
nuclear quadrupole rotational-vibrational structure, from spher- 
ical oscillator to axial rotor to triaxial rotor. Specifically, the 
direct product basis obtained from an optimally chosen set of 
SU(1, 1) p wave functions [17] and the SO(5) D SO(3) spher- 
ical harmonics x Y va LM(y, Q) [4] provides an exceedingly effi- 
cient basis for numerical solution of the Bohr Hamiltonian [5], 
For application to transitional and deformed nuclei, the method 
yields order-of-magnitude reductions in the basis size needed 
for convergence, as compared to diagonalization in a five- 
dimensional oscillator basis. The SU(1, 1) x SO(5) algebraic 
structure of the basis facilitates construction of matrix elements 
for a wide variety of potential and kinetic energy operators. 

In this letter, detailed quantitative predictions for single- 
phonon and multi-phonon excitations in deformed rotor nuclei 



are established by exact numerical diagonalization of the Bohr 
Hamiltonian, making use of newly-calculated SO(5) D SO(3) 
Clebsch-Gordan coefficients [7]. In the past, interpretation of 
rotational phonon states within the Bohr description has largely 
been at a schematic level (e.g., Refs. [18-22]): adiabatic sep- 
aration of the rotational and vibrational degrees of freedom is 
assumed, the (5 and y excitations are taken to be harmonic, and 
phonon selection rules are assumed for electric quadrupole tran- 
sitions. These predictions are then adjusted by spin-dependent 
band mixing [23] with ad hoc mixing parameters. Here, in- 
stead, we explore the actual predictions of the Bohr Hamilto- 
nian. The signatures for the onset of rigid triaxial deformation 
within the Bohr framework are also considered. Preliminary 
results were presented in Ref . [24] . 

The Bohr Hamiltonian [2] is given, in terms of the 
quadrupole deformation variables /3 and y and Euler angles Q, 

by 

h 2 r 1 d , d A 2 1 



H = 

2B 



where A 2 is the angular (y, £2) part of the Laplacian in five di- 
mensions. The essential aspect of the present ACM solutions is 
that the angular degrees of freedom are treated in full, includ- 
ing dynamical y deformation and its coupling to the rotational 
motion. In the context of small-oscillation approximations for 
y (e.g., Ref. [2, 13]), the y dependence of the potential is simply 
taken as oc-y 2 , but for solution of the full problem the y depen- 
dence must be defined more completely. From the symmetry 
properties of a quadrupole-deformed nucleus [10], the poten- 
tial energy V(J3,y) must be periodic in y (with period 120°), 
and it must be symmetric about y = 0° and y = 60°. The 
simplest potential of this form, with minimum at y = 0°, is 
V(y) oc (1 — cos 3y) [3], as shown in Fig. l(a,inset). 

Let us first consider a problem proposed by Iachello [14], 
with Hamiltonian 



H = A 2 + *[(1 - cos 3y) + £cos 2 3y]. 



(2) 
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Figure 1: Level schemes for the angular problem Hamiltonian (2), for (a) f = with^ = 50, (b) £ = 0.5 with^f = 100, and (c) f = 0.8 with^ = 500. The potential is 
shown in the inset, with the ground, quasi-y, and quasi-yy band head energies indicated. (d,e,f) Staggering of level energies within the quasi-y band, for the same 
three calculations, as measured by the energy second difference S(L). 



Only the angular variables (y, D.) are considered, with /3 held 
fixed [25]. The possible forms of the potential V(y) and the re- 
sults of illustrative calculations are shown in Fig. 1. For £ = 0, 
a simple (1 - cos3y) potential is obtained [Fig. 1(a)], approx- 
imately harmonic (locally oc y 1 ) around y = 0°. For £ = 1/2, 
the potential is more softly confining in y, with a quartic min- 
imum (oc y 4 ) at y = 0° [Fig. 1(b)]. This case is termed "criti- 
cal" in Ref. [14]. For £ > 1/2, the potential has a minimum at 
some nonzero value of y, given by cos 3yo = 1 / (2£) [Fig. 1 (c)] . 
The basis functions for the diagonalization consist simply of the 
SO(5) D SO(3) spherical harmonics ^vaLMiy, ^)- The matrix 
elements of physical operators with respect to this basis, both 
for the Hamiltonian and for electromagnetic transitions, can be 
computed directly from the SO(5) D SO(3) Clebsch-Gordan co- 
efficients [4, 7, 26]. Seniority quantum numbers v < 50 amply 
suffice for convergence of the calculations shown. 

The nature of the spectrum obtained depends both on the 
shape of the potential (determined by £) and on the depth of 
the potential (determined by x). The low-lying states form 
quasi-bands which may be roughly identified with the y vi- 
brational excitation (K = 2) and two-phonon y excitations 
(K = and 4). For each calculation in Fig. 1, % lS chosen 
to give E(2y)IE(2\) « 10, appropriate to the well-deformed 
rare earth nuclei. Principal spectroscopic properties consid- 
ered here include the band energies, the detailed level spacings 
within the bands, and the interband electric quadrupole transi- 
tion strengths. 

With the onset of triaxiality (increasing £), the two- 
phonon energy anharmonicities evolve from slightly negative 



(E r y/Ey < 2) for £ = [Fig. 1(a)] to positive (E ry /E y > 2) 
[Fig. 1(c)]. The anharmonicity of the excited K = quasi-band 
rises more rapidly than that of the K = 4 quasi-band. Qualita- 
tively, this is consistent with evolution towards a y-stiff, adia- 
batic triaxial rotor [26, 27], for which the K = 4 quasi-band is 
a triaxial rotational excitation and the K = quasi-band is a y 
vibrational excitation. 

It is essential to observe that the stiffness of the potential 
around its minimum in y determines not only the y vibra- 
tional energy scale but also how well confined the wave func- 
tion is with respect to y. Thus, within the framework of the 
Bohr Hamiltonian, the y band energy [more specifically, the 
energy ratio £(2p/£'(2|), or separation of vibrational and ro- 
tational energy scales] and the y softness of the wave function 
are inextricably linked. From Fig. l(a,inset), it is seen that for 
£'(2p/£'(2[) « 10 the y confinement is weak, and that the range 
of energetically accessible y values increases significantly for 
successive phonon excitations. Confinement is almost nonexis- 
tent at the energy of the two-phonon excitations. 

Consequently, dynamical y deformation plays a major role 
in the calculated structure, as reflected in significant deviations 
from ideal rotational behavior in the spectroscopic predictions. 
Level energies within the y quasi-band [Fig. 1(a)] follow a gen- 
tly y-soft staggering pattern [2(34)(56) . . .]. (The relation be- 
tween the y excitation energy and residual level energy stag- 
gering was noted for transitional nuclei in Ref. [28].) The de- 
viations are even more pronounced for the two-phonon quasi- 
bands. Note especially the near doubling of the rotational con- 
stant, or rotational energy spacing scale, for the two-phonon 
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Figure 3: Interband transition amplitudes, from the y quasi-band to the ground 
state band (top) and from the K = 4 yy quasi-band to the y quasi-band (bottom), 
in Mikhailov form. Values are shown for the calculations of Fig. 1, with f = 
( X = 50) (left), f = 0.5 (x = 100) (middle), and f = 0.8 (x = 500) (right). 
The values shown are for transitions between levels with L < 6, normalized to 
B(E2;2\ -»0+)=l. 



Figure 2: Evolution of spectroscopic properties with y stiffness for the angular 
Hamiltonian (2), for a purely axial potential (f = 0). Quantities shown are 
(a) the energy ratio £ , (4|)/£(2|), (b) excitation energies of low-lying levels, 
normalized to E(2*), and (c) electric quadrupole transition strengths from the 
quasi-y band head to the ground band members, normalized to B{E2\ 2* — > 0|). 



bands relative to the ground state band. 

With increasing £ the level energies within the y band 
progress from y-soft staggering to the pattern associated with 
triaxial rotation [(23)(45) . . .] [12]. This may be seen most 
clearly from plots of the level energy second difference S (L) = 
[[E(L) - E(L - 1)] - [E(L - 1) - E(L - 2)]]/£(2+) [Fig. l(d-f)], 
which has minima at even L for y-soft staggering and at odd L 
for triaxial staggering. As surveyed in Ref. [29], the data for 
most rotational nuclei yield S (L) plots which are either y-soft 
[Fig. 1(d)] or near-constant [Fig. 1(e)]. 

The relation between y softness and spectroscopic properties 
is more systematically and quantitatively examined in Fig. 2. 
For Hamiltonian (2), at fixed the parameter % controls the 
depth and hence y stiffness of the potential. The evolution 
of energy and transition strength prediction with respect to x 
is shown in Fig. 2 for the pure V(y) oc (1 — cos3y) poten- 
tial = 0). Note especially the correlation between the y 
band energy [Fig. 2(b)] and the ground state band energy ra- 
tio E(4*)/E(2*) [Fig. 2(a)], which varies from 2.5 for y-soft 
rotation to 3.33 for rigid axial rotation. This ratio is commonly 
taken as an indicator of rotational adiabaticity. (The quantitative 
details are affected also by the fi degree of freedom.) The evo- 
lution of multi-phonon band energies, approaching harmonicity 
for large %, can also be traced in Fig. 2. The electric quadrupole 
branching ratios, shown for the 2^ state in Fig. 2(c), approach 
the Alaga rule ratios [30] of the adiabatic axial rotor, but only 
slowly, as the y stiffness increases. 

The phenomenological analysis of electromagnetic transition 



Table 1: Effective interband intrinsic matrix elements, as extracted from the 
calculations of Fig. 3. The harmonic axial [27] and Y(5) triaxial [14] estimates 
are included for comparison, along with the experimental values for 162 Dy [31]. 
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strengths in rotational nuclei is founded upon the reduction of 
these strengths to a single intrinsic electromagnetic matrix ele- 
ment and single mixing parameter between each pair of bands, 
according to the Mikhailov mixing formalism [23]. Within 
this framework, all transition amplitudes fall on a straight line 
on an appropriate (Mikhailov) plot, and the intrinsic matrix 
elements and mixing parameter are identified from the slope 
and intercept [2, Sec. 4-4]. The electric quadrupole transi- 
tion strengths for the fully-converged Bohr Hamiltonian calcu- 
lations described above (Fig. 1) are shown in Mikhailov form in 
Fig. 3. Deviations from linearity are significant for transitions 
involving the two-phonon excitations in the £ = calculation 
[Fig. 3(d)], as might be expected from the substantial y soft- 
ness already noted. Otherwise the transition amplitudes at least 
approximately follow an essentially linear pattern. It is there- 
fore meaningful to extract effective intrinsic matrix elements 
and mixing parameters for comparison with experiment. The 
matrix elements for the y — » g transitions and for the yy — » y 
transitions (relative to y — > g) are listed in Table 1 . The adia- 
batic harmonic values [27] and a previous estimate [Y(5)] for 
the "critical" (£ = 0.5) case [14] are shown for comparison. 
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Figure 4: Exact numerical diagonalization of the Bohr Hamiltonian for a well- 
deformed rotor with the full dynamics of the fj, y, and rotational degrees of 
freedom. Confinement in fi is provided by a Davidson potential (fig = 3) with 
y stiffness parameter \ — 5 (see text). (Left) Level energies, normalized to 
£(2|), with an illustration of the Davidson potential (inset). The dashed lines 
indicate the asymptotic forms fi and /? of the potential, at small and large 
fi, respectively. (Right) Transition amplitudes: (b) y — » g, (c) fi — > g, and 
(d) fi — > y. The values shown are for transitions between levels with L < 6, 
normalized to B(E2; 2\ -> Op = 1. 



The nucleus 162 Dy has recently been the subject of detailed 
spectroscopic study by Aprahamian et al. [31], yielding exten- 
sive sets of experimental values for interband electric transition 
strengths. The experimental intrinsic matrix elements are given 
in Table 1 . The y band and candidate two-phonon y band en- 
ergies for 162 Dy closely match those of the £ = (x = 50) cal- 
culation of Fig. 1(a). The measured y — » g strength is about 
a factor of two lower than the Bohr Hamiltonian predictions, 
and the K = 4 yy — > y strength is yet a factor of two lower 
again, i.e., even taken relative to this already reduced y — > g 
strength. Moreover, the observed level energies in 162 Dy [31] 
conform much more closely to adiabatic rotational L(L +1) 
spacings than expected from the Bohr Hamiltonian calculation 
[Fig. 1(a)]. Although we do not include a detailed quantitative 
analysis here, for the y band staggering, compare Fig. 3(d) of 
Ref. [29] with the present Fig. 1(d). The excited bands in 162 Dy 
are actually observed to exhibit decreased rotational constants, 
relative to the ground state band, in stark contrast to the present 
calculation (a similar discrepancy is noted [20] in comparison 
of interacting boson model [32] predictions with data). A more 
detailed comparison requires exploration of the interaction of 
the /3 and y degrees of freedom. 

Let us briefly examine the results of a calculation involving 
the full dynamics of the Bohr Hamiltonian (1), including the 
j3 degree of freedom. An essentially unlimited variety of com- 
bined (3 and y dependences {e.g. , ff 1 cos" 3y) may be used in the 
ACM potential [26]. For illustration, we consider a Davidson 
potential for (3 [Fig. 4(a,inset)], together with the simplest ax- 



ial y confinement from above, so V(J3, y) oc (J3 Q /fi - /3/fio) 2 + 
X(l - cos3y). Without loss of generality, we may take the 
proportionality constant and the inertial constant h 2 /(2B) to be 
unity, if only energy and E2 strength ratios are to be consid- 
ered [33]. This leaves freedom only in the stiffness (fio) and 
y stiffness (x). A fully converged ACM calculation, for /3o-3 
and x = 5, is shown in Fig. 4(a), with quasi-band assignments 
indicated. Fewer than five SU(1, 1) basis functions are required 
for convergence. Electric quadrupole transitions connecting the 
ground,/?, and y bands are shown in Fig. 4(b-d). The zero in the 
Mikhailov plot for the -» y transitions [Fig. 4(d)] reflects the 
zero expected [at £2(^2 + 1) - ^1(^1 + 1) =4] if simple phonon 
selection rules hold on the intrinsic matrix elements (see Fig. 12 
of Ref. [34]). 

The possibility of exact diagonalization of the Bohr Hamil- 
tonian for essentially arbitrary f3 and y stiffness, by means 
of the algebraic collective model, opens the door for direct 
comparison of the Bohr Hamiltonian predictions with experi- 
ment throughout the range of possible dynamics for the nuclear 
quadrupole degree of freedom. At a phenomenological level, 
this permits meaningful tests of the Bohr Hamiltonian for gen- 
eral rotor-vibrator nuclei. More fundamentally, in the context 
of nuclear many-body theory, it is essential to know the limita- 
tions and necessary modifications to the Bohr framework, since 
microscopic descriptions of nuclear collectivity rely upon a re- 
duction of the many-body problem to one involving effective 
collective degrees of freedom, typically those of the Bohr de- 
scription [35, 36] or its symplectic generalization [37]. 

Notably, these preliminary results suggest that the Bohr de- 
scription quantitatively overpredicts and even qualitatively mis- 
predicts the nature of deviations from adiabatic rotational struc- 
ture, and that these deviations result from dynamical y softness 
inherently incorporated into the description. The Bohr descrip- 
tion also strongly overpredicts y phonon interband transition 
strengths (although this could already be anticipated from, e.g., 
rotation-vibration model estimates). A basic question is there- 
fore the extent to which the discrepancies lie in the details of 
the description {e.g., restriction of the kinetic energy operator 
to quadratic order in the momenta) or in a fundamental failure 
of the collective quadrupole degree of freedom to account for 
the observed phenomena. 

Valuable discussions with N. V. Zamfir, D. J. Rowe, 
F. Iachello, S. Frauendorf, and A. Aprahamian are gratefully 
acknowledged. This work was supported by the US DOE un- 
der grant DE-FG02-95ER-40934. 
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